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Abstract This paper is devoted to the analysis of non-negative solutions for a de- 
generate parabolic-elliptic Patlak-Keller-Segel system with critical nonlinear diffusion in a 
bounded domain with homogeneous Neumann boundary conditions. Our aim is to prove 
the existence of a global weak solution under a smallness condition on the mass of the 
initial data, there by completing previous results on finite blow-up for large masses. Un- 
der some higher regularity condition on solutions, the uniqueness of solutions is proved by 
using a classical duality technique. 
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1 Introduction 

Chemotaxis is the movement of biological organisms oriented towards the gradient of some 
substance, called the chemoattractant. The Patlak-Keller-Segel (PKS) model (see [13], [12] 
and [17]) has been introduced in order to explain chemotaxis cell aggregation by means of 
a coupled system of two equations: a drift-diffusion type equation for the cell density u, 
and a reaction diffusion equation for the chemoattractant concentration cp. It reads 



(PKS) < 



d t u = div(Vu m - u ■ Vip) x € ft, t > 0, 

-Aip = u- < u > x € ft, t > 0, 

< <p(t) > = t > 0, (1) 

d„u = d u ip = x G 9ft, t > 0, 

u(0,x) = Uq(x) x G ft, 



where ft C M. N is an open bounded domain, v the outward unit normal vector to the 
boundary 9ft and m > 1. An important parameter in this model is the total mass M of 
cells, which is formally conserved through the evolution: 

M =< u >= — — / u(t, x) dx = — ■ I uq(x) dx. (2) 

Several studies have revealed that the dynamics of (1) depend sensitively on the parame- 
ters N, m and M. More precisely, if N = 2 and m = 1, it is well-known that the solutions 
of (1) may blow up in finite time if M is sufficiently large (see [17, 16]) while solutions are 
global in time for M sufficiently small [17], see also the survey articles [4, 10]. 



The situation is very different when m = 1 and N ^ 2. In fact, if N = 1, there is 
global existence of solutions of (1) whatever the value of the mass of initial data M, see 
[8] and the references therein. If N > 3, for all M > 0, there are initial data uq with mass 
M for which the corresponding solutions of (_1) explode in finite time (see [16]). Thus, 



in dimension N > 3 and m = 1, the threshold phenomenon does not take place as in 
dimension 2, but we expect the same phenomenon when N > 3 and m is equal to the 
critical value m = m c = 2 ^ N r^~ ■ More precisely, we consider a more general version of (V) 
where the first equation of (1) is replaced by 

dtu = div((f)(u) Vu — u Vtp), t > 0, x G fi, 

and the diffusitivity <fi is a positive function in C 1 ([0,oo[) which does not grow to fast 
at infinity. In [8], the authors proved that there is a critical exponent such that, if the 
diffusion has a faster growth than the one given by this exponent, solutions to (1) (with 
4>(u) instead of mii m_1 ) exist globally and are uniformly bounded, see also [6, 14] for 
N = 2. More precisely, the main results in [8] read as follows: 

• If <j){u) > c(l + u) p for all u > and some c > and p > 1 — jj then all solutions of 
(1) are global and bounded. 

• If 4>{u) < c(l + u) p for all u > and some c > and p < 1 — A then there exist 
initial data tto such that 

lim | t) | |cso = for some finite T > 0. 

t— >T 

Except for N = 2, the critical case m = 2 ^ j ^~ 1 ^ is not covered by the analysis of [8]. Re- 
cently, Cieslak and Laurengot in [7] show that if <j){u) < c(l + u) l ~^ and > 3, there are 
solutions of (1) blowing up in finite time when M exceeds an explicit threshold. In order 
to prove that, when N > 3 and m = 2 ^ — , we have a threshold phenomenon similar to 
dimension N = 2 with m = 1, it remains to show that solutions of (1) are global when M 
is small enough. The goal of this paper is to show that this is indeed true, see Theorem 
2.2 below. 

By combining Theorem 2.2 with the blow-up result obtained in [7], we conclude that, 
for N > 3 and m = ja7 , there exists < M\ < M<± < oo such that the solutions of (1) 
are global if the mass M of the initial data uq is in [0, Mi), and may explode in finite time 
if M > M%. An important open question is whether Mi = M2 when fi is a ball in R N and 
uq is a radially symmetric function. Notice that, in the radial case, this result is true when 
N = 2 and m = 1, and the threshold value of the mass for blow-up is Mi = M% = 8tt, 
see [6., I6j 15., 18]. Again, for A^ = 2 and m = 1, but for regular, connected and bounded 
domain, it has been shown that Mi = 4-7T = (see [15 . 16] and the references therein). 

Such a result does not seem to be known for A^ > 3 and m = 2(,N N ^ . 

Still, in the whole space 17 = when the equation for ip in (1) is replaced by the 
Poisson equation ip = En * u, with En being the Poisson kernel, it has been shown in 
% 5, 2, 20, 21, 3] that: 

• When A^ > 3 and 1 < m < 2 — j^, this modified version of (1) has a global weak 
solution if M = ||ito||i is sufficiently small, while finite time blow-up occurs for some 
initial data with sufficiently large mass. 

• When N > 2 and m > 2 — , this modified version of (1) has a global weak solution 
whatever the value of M. 
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• When N > 2 and m = 2 — ■£>, there is a threshold mass M c > such that solutions 
to this modified version of (JL) exist globally if M = ||ito||i < M c , and might blow up 
in finite time if M > M c . 

From now on, we assume that 

2(N-1) 

N>3 and m = — -. 

N 

2 Main Theorem 

Throughout this paper , we deal with weak solutions of (1). Our definition of weak 
solutions now reads: 

Definition 2.1. Let T £ (0;oo]. A pair (u,(p) of functions u : £1 x [0, T) — > [0, oo), 
ip : £1 x [0, T) — > R is called a weak solution of (1) in £1 x [0,T) i/ 

• u £ L°°((0, T); L°°($7)); n m G L 2 ((0, T); (0)) and < n >= M. 

• if e L 2 ((0,T);H x (p)) and <tp>=0. 

• (u,ip) satisfies the equation in the sense of distributions ; i.e, 

- [ [ (Vu m ■ VV> - uVip ■ VV> - u d t ip) dxdt= [ u Q (x) ip(Q, x) dx, 
Jo Jn Jn 

[ [ V93 • dxdt = f [ (u-M) dxdt, 
Jo Jn Jo Jn 



for any continuously differentiable function ip £ C 1 ([0, T] x Q) with ip(T) = and T > 0. 
For ip G H 1 (Q) satisfying < cp >= 0, we denote by C s the Sobolev constant where 

||V^|| 2 >C S |H| 2 ., where 2* = (3) 

The main theorem gives the existence and uniqueness of a time global weak solution to 
(1) which corresponds to a degenerate version of the "Nagai model" for the semi-linear 
Keller-Segel system, when uq £ L°°(f2) and the initial data is assumed to be small. 



JV 
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M*:=\ - 2 5,„, a 1 , (4) 



Theorem 2.2. Define 

M* := I 

(m — 1) \Q\n 
where C s is the Sobolev constant in (3). 

Assume that uq is nonnegative function in L°°({1), which satisfies 

|K||i<M*. (5) 

Then the equation (1) has a global weak solution (u, <p) in the sense of Definition 2. 1. 
Moreover, if we assume that 

^GL oo ((0,T);VF 2 ' oo (^)) (6) 
for all T > then this solution is unique. 
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In order to prove the previous theorem, we introduce the following approximated equa- 
tions 



(KS) S { 



d t u s = div (V(u s + 5) m - u s X7ip s ) x G O, t > 0, 

— Aips = us— < us > xG$7,t>0, 

d v us = d v ips = x G dU,t > 0, 

ug(0,x) = uq(x) x G f2, 



where 5 G (0, 1), and we show that under a smallness condition on the mass of initial data, 
the Liapunov function 

L s (u,<p) = / (6,5 (it) + \\Vips\ 2 - u s <ps) dx, 
Jn 1 

yields the L m bound of ug(t) independent of 5. Then using Gagliardo-Nirenberg and 
Poincare inequalities, we obtain for p > m, the L p bound for ug{t) independent of 5. 
As a consequence of Sobolev embedding theorem, we improve the regularity of ips- And 
thus, under the same assumptions on the initial data, Moser's iteration technique yields 
the uniform bound of us- Then, thanks to the local well-posedness result [8, Theorem 
3.1] we obtain the existence of a global solution of (KS)s- The existence of solutions 
stated in Theorem 2.2 is then proved using a compactness method; for that purpose we 
show an additional estimate on dtu™ which, together with the already derived estimates, 
guarantees the compactness in space and time of the family (1*5)5^(0,1)- Finally, in the 
presence of nonlinear diffusion and under some additional regularity assumption on (ps, we 
prove the uniqueness using a classical duality technique. 



3 Approximated Equations 



The first equation of (JL) is a quasilinear parabolic equation of degenerate type. Therefore, 
we cannot expect the system (1) to have a classical solution at the point where u vanishes. 
In order to prove Theorem 2.2., we use a compactness method and introduce the following 
approximated equations of (KS): 



(KS)s { 



d t us 
-Aips 
d v us = d u ips 
ug(0,x) 



div(V(ii5 + <5) r 
us- < u s > 


U (x) 



usVips) 



x G Q,t > 0, 

ie!l,f>o, 

x G dQ,t > 0, 
x G fi, 



(7) 



where 5 G (0, 1). 

The main purpose of this section is to construct the time global strong solution of (7). 



3.1 Existence of global strong solution of (KS)g 

Theorem 3.1. For 5 G (0,1) and T > 0, we consider an initial condition uq G L°°(£l), 
uq > and such that ||ito||i < M* where M* is defined in (4). Then (KS)g has a global 
strong solution (ug,(pg) which is bounded in L°°((0,T) x f2) for all T > uniformly with 
respect to 5 G (0, 1). 

The starting point of the proof of Theorem 3.1 is the following local well-posedness 
result [8, Theorem 1.3]: 
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Lemma 3.2. Let the same assumptions as that in Theorem 3.1 hold. There exists a maxi- 
mal existence time T^ ax . € (0, oo] and a unique solution (us, ps) of (KS)s in [0, T^ ax ) x O. 
Moreover, 

if T^ax < oo then lim Hn^t, -)lloo = oo. 

t^T s 

" ' max 

In addition < u$(t) > = < no >= M for all t £ [0,T^ ax ). 

To prove Theorem 3.1 we need to prove some lemmas which control L m norm, LP norm 
and L°° norm of the solution us of (7). 



3.2 L p -estimates, 1 < p < oo. 

Our goal is to show that if ||no||i is small enough then all solutions are global in time and 
uniformly bounded. 

Let us first prove the L m bound for us- 

Lemma 3.3. Let the same assumptions as that in Theorem 3.1 hold and (us,ps) be the 
nonnegative maximal solution of (KS)s- Then, us satisfies the following estimate 

\\us(t)\\ m < C , for all t G [0,T max ) 

and \ \us(t)\\i = ||no||i where Cq is a constant independent ofT^ ax and 5. 

Proof. In this proof, the solution to equation (7) should be denoted by (us,ps) but for 
simplicity we drop the index. 
Let us define the functional Ls by 

L 5 (u, p) = [ (b s (u) + ^|Vi^| 2 -u ip) dx, 
Jn z 

where 



h Ji a 

such that 6(5 (1) = b' s (l) = and b(u) > 0. According to [U] it is a Liapunov functional 
for (KS) S . Indeed, 

—Ls(u(t),(p(t)) = j b' s (u) dtu dx — ( Ap dtp dx — ( dtu p dx — j u dtp dx 
dt Jr, In Jo Jo 



dtu (b's(u) — ip) dx — j (Ap + u) dtp dx 
Jn 

/ div (m (u + <5) m_1 Vu-u Vp) (b' s (u) - ip) dx- I < u(t) > d t p dx 
Jn Jn 

- [ (m(u + 5) 771 - 1 Vu-u Vp) (bg(u) Vu — Vp) dx - M ^- [ 
Jn dt J n 

u (b's(u) V-u - Vp) 2 dx 



p dx 



< 0. 

Then, we can conclude that for all t £ [0, T„ iax ) we have Lg(u(t), p(t)) < Ls(uq, po). Using 
Sobolev inequality (3), Holder inequality, and Young inequality we obtain 



1 C~ 2 

U p dx < \\p\\2* |M| 2JV < C~ 1 \\Vp\\2 \\u\\ 2JV < -\\Vp\\2 H fr~\\ V ~->.\ • 



|2 _|_ n .. 2 

JV+2 ~~ S M ' 1 1 - 11 jv+2 _ 2 " r ' 2 JV+2 
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Since < m i an d using interpolation inequality we get, 

1_ N-l 1 

\\u\\^<\\u\\i \\u\\ m N <Mn \n\77 \\ u 

Then, 



2 

m ■ 

N+2 



Substituting this into the Liapunov functional, we find: 



J^u<pdx< ^\\V<p\\ 2 2 + ^- M# ||u| lm . 



L S {u,<p) > f(bs(u) + ^\V^\ 2 )dx-h\Vip\\ 2 2 -^j-M^\n\^\\u 

> f bs(u) dx - M^f [|u||™. 
Jn 2 



We next observe that: 

J — dsdz - m J J sm ~ 2 dsdz 



Then: 



u m u m 

> - -u + l> -u. 

m — 1 m — \ m — 1 m — 1 



m—\ I m—\ 

1 C7 2 , 2 ,„, 2 \ ,, ,, m m 



m-1 2 / m-1 



Miv \\ u \\™ -M\Q\. 



Let us define ojm by 



x-.w — - - ^- .WT = (M* -MW). 

in — 1 Z z u s 

Since M = ||wo||i < then ojm is positive. Finally we get, 

L s ( Uo ,ip ) + -^M\n\>L s (u(t),ip(t)) + -^M\n\>u; M \\u(t)\\™ forte [0,T^ o J. 

In addition, we can see that Lg(uo, (fo) < C where C is independent of 8 € (0, 1). In fact, 

Ls(u ,(p ) = [ (b s (u ) + hx7(p \ 2 - u po) dx, 
Jn 1 

and, since (5 + s) m ~ l < 5 m ~ l + s m ' 1 < 1 + s™" 1 we obtain 

p f QS + g)™- 1 p f l + s"^ 1 

A A « A Ji s 

m ( Uq 1 \ 

< m{u In u - u + 1) H u + 1 . 

m — 1 \ m J 

Using Young inequality we get 

t i \ I 1 1 2 i^i Ikollm mlfil l.._ ll2 1,, ,,2 1,, ,,9 

-M«o,<A)) < m 1 1 zio 1 1 2 + m l"l H — 7 H — r + dlv^olb + dpolb + dl^olb- 

TTl X 777) J. Zi Zi Z 

since uo € L°°(Q) and <^o € H l (Q) we get Ls(uo,ifo) < C where C is independent of <5 
and the proof of the lemma is complete. □ 
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Thanks to Lemma 3.3, let us now show that for all p > m the L p bound for ug. 

Lemma 3.4. Let the same assumptions as that in Theorem 3.1 hold. Then for all T > 
and all p G (l,oo) there exists C(p,T) independent on 5 such that, for all t G [0,T^ ax ) n 
[0,T], the solution (us,<ps) to (KS)g satisfies 

\\u s (t)\\ p <C(p,T), (8) 

and 

ft r 

(5 + n 5 ) m ^ 1 u p 5 ' 2 \Wu s \ 2 dxds < C(p,T). (9) 



'o Jn 

To prove the previous lemma we need the following preliminary result [20]. 
Lemma 3.5. Consider < qi < q2 < 2*. There is C\ depending only on N such that 

\\u\\ q2 < C{ [Mll^n) \\u\\\-\ for u G H\n), (10) 

with 

2Af fe-«) el0 ,l]. 



« 2 [(JV + 2) gi +2JV(l- gi )] 
Proof. For u G H 1 ^) we have by Sobolev inequality 

1Mb* < Civ||u[|_H-i. (11) 
By interpolation inequality we have for < qi < q2 < 2* 

IHIw < ll«[|2. (12) 

where ^ = d ^ N 2 ^ + Hence, substitute (11) into (12) and the lemma is proved. □ 

Now, we recall the following generalized Poincare inequality. 
Lemma 3.6. For u G H l (£l) we have for < q\ < 1 the following inequality 

\H\ 2 m <C 2 {qi) (llVulH + IMgj, 

where C2 depends only on f2 and q\ . 

Now using the last two lemmas, let us prove Lemma 3.4 . 

Proof. In this proof, the solution to equation (7) should be denoted by (us,(ps) but for 
simplicity we drop the index. 

We choose p > 1, K > and we multiply the first equation in (7) by {u — K) p + l and 
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integrate by parts using the boundary conditions for u and tp to see that 
^\\( U -K)+\\P = -m{p-l) j^5 + u) m - 1 {u-Kf- 2 \Vu\ 2 dx 



+ (p - 1) / u V(/? (u - X)+ 2 • V« ds 

= -m(p -1) / + K + i^)" 1 - 1 (u - Kf- 2 |Vu| 2 dx 

+ - 1) / (u - K + K) V<p ■ (u - K)+~ 2 Vu dx 
Jn 

< -m(p - 1) f (u + 5- K) 171 - 1 (u - K) p + 2 |Vu| 2 dx 

Jn 

+ (p - 1) / (u - V( /? • Vn dx + (p - f Vip (u- K) p - 2 ■ Vu 

Jn Jn 

< -m(p - 1) f (u + 5- K) m - X (u - Kf- 2 |Vu| 2 dx 

Jn 

- — 1 / (« - Kf + Aipdx-K ( (u- Kf- 1 Aip dx 
Jn Jn 



dx 



V Jn 

< —m 



i(p - 1) f (u + 5- K)" 1 - 1 [u - Kf- 2 \Vu\ 2 dx + (J), 
Jn 

where, thanks to the second equation in (7), 

(/) = ^— !- f (u - K)l (u - M) dx + K [ (u - Kf' 1 (u - M) dx 
P Jn Jn 

= - K) + \\lH + P —^(K - M)\\(u - K) + \\l 

+ K\\(u - K)+\\' + K(K-M)\\(u- K) + \f p Z\ 

< K 2 \\{u - K) + \\IZ\ + 2K \\(u - K) + \\? + \\(u- K)+\%1\. 

Since for a > and b > we have a p ~ l b < a p+1 + 6^ and a p b < a p+1 + then, 

(/) < 2,\\{u-K) + \\ p p X\ + {2K) p ^+K p+ \ (13) 

and we get 

^-\\{u-K) + \\ p < -m(p - 1) / (u + <5 - K)™- 1 (u - K) p ~ 2 \Vu\ 2 dx 
dt ' Jn 

+ 3p\\(u-K) + \\ p p Xl + C p K p+ \ (14) 

for all ie [0,3^). 

The term ||(u — K) + 11?+} can be estimated with the help of Lemma 3.5 and Lemma 3.6. 
Assuming now that p > 2 we remark that < < 1 and 1 < Jx^zi = jfz% 1+ n p < 

2N then thanks to Lemma 3.5 and Lemma 3.6 we obtain 



N-2> 



p+m — 1 

2(p+l) 

/ p+m — 1 2(p+l) a p-\-m — l 

{u-K) + * < C(p) ||V(u-K) + 2 HI+- 1 2 II' 1 '" J - 

p+m — 1 \ 



2 



p+m— 1 
p+m — 1 

+ h(u-k) + 2 irr 1 , as) 

p+m— 1 J 
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where 

= ^f±e(o,l). (16) 

Since 

p+m-l 2(p+i) /" 

ll(n-^) + 2 H3£ = / (n-ETf + +1 dx=||(n-Er) + ||^, (17) 

p+m — 1 i/ Q 

p-f-m — 1 

*£±lL n a\ / f \ (P+l)(l-f) 2 

\\{u-K) + * \\^TT } = ( (u-K) + dx) =\\(u-K) + \\», (18) 

and by Lemma 3.3 

1 f Tsm-l „, j„ ^ \H\m / Co 

, /v ^ m ~ l Ju>K 



u-K) + \\ 1 = (u-K)dx<- =-/ K" 1 - 1 udx<^P^r<-P—, (19) 



we substitute (17), (18) and (19) into (15) and obtain 

\\{u-K)+\ \0 < C 3 (p) 1 1 1 V(u - K)f^ ||1 K~~^r^- + ^ (m-i)<P+U | (20) 
We may choose = large enough such that 



-2(m-l) 



3pC 3 (p)if* " < 



2 p (p — 1) m 
[m + p — l) 2 ' 



Hence 

so that 
As 



^||(u-^) + |!£<c(p)^+\ 



1 1 («(*)-#»)+! IS <C(P) t+||«o||5, for ie [0,7^' 



5 ) 

max / 



[ \u\ p dx < [ (2 K.f- 1 \u\ dx + [ \u - K* + K*\ p dx 

Jn Ju<2K t Ju>2K* 

< {2K tf y- 1 M+l {2 \u- K*\) p dx, 

Ju>2K t 

< (2KJP- 1 M + 2P \\{ U -K,)+\\*, 
the previous inequality warrants that 

\\u(t)\\ p <C(p,T), t€ [0,T max )n[0,T], (21) 
where C(p, T) is a constant independent of 5. 

We next take K = in (14), integrate with respect to time and use (8) to obtain 
(9). □ 

Thanks to Lemma 3.4, we can improve the regularity of (ps- 

Lemma 3.7. Let the same assumptions as that in Theorem 3.1 hold, the solution ip§ 
satisfies 

\\VMt)\\oc < L(T), te[o,r max )n[o,T] 
where T > and L is a positive constant independent of 5. 



Proof. Using standard elliptic regularity estimates for cp§, we infer from Lemma 3.4 that 
given T > 0, and p G (1, oo), there is C(p, T) such that 

|M*)llv^ < C(p) \\u 5 {t)\\ p < C(p,T), for t G [0,T max ) n [0,T]. 

Lemma 3.7 then readily follows from Sobolev embedding theorem upon choosing p > 
N. □ 

Lemma 3.8. Let N > 3, r > 4, u G L^(U), and u' + ™ 1 G ff^fi). T/ien if Zio/ds i/iaf 



||«||r < Cf m " 1 ||« 2— HSi"- 1 (22) 

4 

wif/t 



3iV(r+m-l) , . 

G(0,1). (23) 



(3JV + 2)r + 4iV(m - 1) 
Proof. For r > 4, we can see that 

i 

r+m-l 2r N 

(U 2 J r + m-l rfj; 


and 

r 2r 2iV 
< 1 < < 2 < 



2(r + m-l) r + m-l N - 2' 

By Lemma 3.5, 



. T / „/) . . r + m — 1 /) r+m — 1 -i /i \ r + m 1 

||U|| r = \\U 2 " £T < Cf U - 5— ^i fm tlT- M r 

r+m-l V 2(r+m-l) 

and 

2/V C — - 1 

Q Z,JV Vr+m-l 2 (r+m-l)' 1 



2r 



1 ( 2 ^(! " 217+^17) + + 2) 2-TT+^T)) 
6(0,1). 



r+m 

3iV (r + m-l) 



(3iV + 2) r + AN (m - 1 
In addition, we have 

2(r+m-l) 



r+m — 1 

1 2(r + m-l) 



r+m — 1 

// - | — ( I \u\ 5 2(r + m-l) _ ~ 



and we obtain (22). □ 

We are now in a position to prove the uniform L°°(0,) bound for us- 

Lemma 3.9. Let the same assumptions as that in Theorem .3.1 hold, and (ug,(ps) be the 
nonnegative maximal solution of (7). For all T > 0, there is Coo(T) such that 

\\u 5 (t)\\oo < Coo(r), for all t G [0,T* ax ) n [0,T], 

where Coo(T) is a positive constant independent on 5. 
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Proof. In this proof we omit the index 5, and we employ Moser's iteration technique 
developed in [1, 21] to show the uniform norm bound for u. 

We multiply the first equation in (7) by u r ~ l , where r > 4, and integrate it over £1. Then, 
we have 



= - [ (V(u + 5) m - u Vif) ■ W -1 dx 
Jn 



dt r 



-m(r - 1) / (u + <5) m_1 u r ~ 2 \Vu\ 2 dx + (r - 1) / u 7 " 1 Vip ■ Vu dx 
Jn Jn 



< —m(r — 1) / u 



( u m+r - 3 \S/u\ 2 dx + (r - 1) ( u r_1 Vip ■ Vu dx. 
Jn Jn 



By Young's inequality and Lemma 3.7, 

Id,, „_ -Amir -I) f .„ r+ m -i ,, , 2(r - 1) NVdU , .,_ , 

— r u I < - i > I |Vn^^| 2 dx + -A- ; M y n°o / ^ 2 |Vn^^| dx 

r dt (r + m — l) A Jq (r + m — 1) 



-4m(r 


-1) 


(r + m - 


-I) 2 


—2m(r 


-1) 


(r + m - 


-I) 2 


—2m(r 


-1) 


(r + m - 


-If 



n 



r— m+1 

u 

^ 7^ Tjl|Vu^"||^ + C(r) r / u r ~ m+1 dx. 



n 

Using Holder and Young inequalities and Lemma 3J3 we obtain 

r dt (r + m — l) z 

< (r "^ l|Vn^||| + C- + r 2 |H|;, (24) 

(r + m-1) 2 

r-1 - 

where we have used that r r - m < r for r > 4. 
By Lemma 3.8, we have for r > 4 



| u ||r < C rlL-i || u |[Ki-») ||uT||^-\ (25) 



where 

3iV(r + m-l) 



(3JV + 2)r + 4iV(m - 1) 
Therefore, Young inequality and (25) yield that 



< 1. 



2r 2 \\u\\ r r < 2r 2 \\u\\f-°\\u^\\^ 1 

4 

Or m (r — 1) r + m — 1 r+m-i 2 

~ r + m-1 (r + m-1) 2 r C 2 (l) " 

9r 

r + m — 1 — Or I 0(r + m — l)r\ r(i-fi)+m-i 

+ ; 1 ^V 1 ) 7 TV 

r + m — 1 \ m(r — 1) 

4 

where C2(l) is the Poincare constant defined in Lemma 3.6. Then we obtain 

m (r — 1) r + m-1 



<-» 2 II I it- ^ " t V ~~ x ) II r + m-i n 



C 2 (l) (r + m-1) 



(r+m-1) 2(r + m-l)+9r 



+ c r(l-6) + m -l 2 r(l-fl) + m-l r r(l-e) + m-l 

4 



(l-8)r(r + m-l) 
I r(l-fl) + m-l 
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Now, since N > 2, which gives 4N > 3N + 2, we find the following upper bound for 9 

< ™ (26) 

~ 3N + 2 v ; 

In addition, 

9 r 9 1 3N 

r(l-0) + m-l ~ 1-9 1-0-2' v; 

r + m — 1 r + m — 1 1 3iV + 2 . 

< - — < < — , (28) 



r(l - 0) + m - 1 ~ (l-9)(r + m-l) ~ 1 
and 

^±^!±^<^<9iV + 4. (29) 
r{l-9) + m-l ~ 1-9 ~ v ' 

As C\ > 1 and r > 1, we get 

/ i \ (l-9)r(r + m-l) 

9r 2 ll7/ll r < m[r-ij r±rpl ||2 r 97V+4 m || r(i-e )+m -i /nM 

2r ||u|| r < (r + m _ 1)2 l|n Ibi + G r |M|r ■ (30) 

Using Lemma 3.6 we have 

r + m— 1 ,,o , > /.. r + m— 1 ,,o . r+m — 1 ,,o\ 

||«^Hhi <0 2 (1) (||V«— 3— + . (31) 



Using Holder inequality, Young inequality and Lemma 3.3, we get 



r + m — 1 



r+m — 3 r — m+1 r+m — 3 

|2 it nr+ra— 1 ^ IL.Il' r-1 Ik. II j — 1 il.m r ~ 1 IL. II r_1 



If 2 III = IMI vI;"-! < IMk IMIi <Nir n«om 

2 

then, 

m (r — 1) r+m-i o r-i r + m - 3 .. ,,_ 

1 1 IX 2 |[f < (r— l)r+m-3 \\ U \\r 



r-1 



(r + m — l) 2 r — 1 

2-m / m n n^^^A ~ 

+ — r i« K i 

r — 1 \ (r + m — 1) 



_r_-l 

™ r+m— 1 \ 2- 

2 I I . I I r- i / 171 II II r-1 x 



< r IMIr + 7 : TT2 IKUl 

\{r + m — l) z 

< r 2 \\u\\ r r + C r A . (32) 

Now substituting (32) and (31) into (30) we get 



I 1 \ , (l-fl)r(r + m-l) 

m V ~ 1 ) i / ||v7., '' + '?~ 1 l|2 _l IU, r+ ? -1 ||2^ _l r< ^9A+4 ||„,|| r(l-0) + m-l 

(r + m — l) 2 



2r^||n||; < ^ ||Vw i ^||^ + ||« i ^||J + C r 9JV+4 \\u\\ L 



4 

(l-fl)r(r+m-l) 



/ -i \ (1 — 6)r(r+m — ] 

< m ( r ~ IL || Vu ^|| 2 + r 2 |Hi; + ci + Cr 9W+4 |H|7™^ , 

(r + m — l) z . 4 

hence 

I 1 \ r (l-fl)(r+m-l) 

r 2 | M |r< m(r-l) ^ 2 ^+4 | | u | | r(l-«)+m-l 

11 llr ~ (r + m - l) 2 11 4 " M 4 

We apply Young inequality again to the last term of the above inequality. It is easy to see 
that 

2 <1 e< (1-9) (r + m-1) _ (1 - 9)r + (1 - 9){m - 1) ^ 



3iV + 2~ ~ r(l-0) + m-l r(l-0) + m-l 
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so that 

AMI < ^^-'h lVu^Wl + CI + 1 + (C r ^) m+1 |Hr P) (33) 
(r + m — l) z 4 

for any r £ [4, oo). 

Substituting (33) into (24) we end up with 

^-\\u\\ r r <r Cl + r + r (C r m+i ) 3N+1 \\u\\ r r < CI + Cr a \\u\\ r r , (34) 

Qt 4 4 

for any r G [4, oo), where a = (9N + 4)(3iV + 1) + 1. After integrating (34) from to t, 
we obtain the L r estimate for u as follows: 

sup \\u{t)\\ r r < \ \u \\ r r + T Cl + C r a T sup \\u{t)\\l. (35) 

0<t<T 0<t<T 4 

Since 

7 1 1 

\\uo\\r < l!'Wo|!oo r Ikolli < C 6 , 

then 

sup \\u(t)\\ r r <C 7 (T) r a max ice, sup , (36) 

0<t<T I 0<t<T 4 J 

and we obtain for r > 4 

sup ||u(t)|| r < C 7 (T)~ r~ max < Cq, sup ||ti(t)||i: > . (37) 

0<t<T I 0<t<T 4 J 

We are now in a position to derive the claimed L°° estimate. To this end, we set 

a p : = max <Cq, sup ||w(i)||4P 
I 0<t<T 

for p > 0. Then we take r = 4? with p > in (37) which reads 

ap < C 7 (T)^ max^C 6 , sup ||u(t)|| 4 p-i > , 

I o<t<T J 

< 4# C 7 (r)4T ap _ 1 
since p < 2 P for p > 1. Arguing by induction we conclude that 

a p <^n=^ k C 7 {T)n=i±- k ao . 
Then by using Lemma 3.3 we get 

sup ||«(t)||4P < 4 a C 7 (T) a < C 8 {T). 

0<t<T 

Consequently, by letting p tend to oo, we see that u G L°°((0,T) x $7) and 

sup Halloo <C 8 (T). (38) 

0<t<T 

Since the right hand side is independent of 5, we have proved the lemma. □ 
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Lemma 3.10. Let the same assumptions as that in Theorem 3.1 hold, and (ug,(ps) be the 
solution to (7). Then for all T > there is Cg(T) such that the solution us satisfies the 
following derivation estimate 



T 



ll^' u ™ll(W /1 ' JV + 1 )' dt < Cg(T). 

Proof. Consider ijj G W 1>JV+1 (n) and t G (0,T), we have 



./ 


m 


u 




l 




m 


/ 


V 




Jn 





m— 1 / 



,m—l 



= m 

< m 

+ 1 

< m 



Jn 



f- 1 \Vuf\ \V^\+uf |X7-0| |V W | 
m(m - 1) uf 1 ^ \Vu s \ 2 + |^| (m - |Vu 5 | |V^|] dx 



l^llST 1 IIVuHh l|vvi| 2 + ||vv|| 2 ||«i||Sllv w ||oo|n|* 



+ IM 



4m(m — 1) 
(2m- l) 2 



|Vn A 



+ 



2 HLJ:||Vu?ll2 ||v w | 
m 



Using Lemma 3.8, Lemma 3.9, and the embedding of W 1 ' in L°°(Q), we end up 

with 



< 



d t uT(t),i>>\<c(T)(\\Vu s (tr\\ 2 + \\Vu™ Hml + ii \m\w^, 

and a duality argument gives 

\\d t uT(t)\\(w^y < C{T) (\\vuf(t)\\ 2 + \\vu^{t)\\l + 1 



Integrating the above inequality over (0, T) and using Lemma 3.4 with p = 2 and p = m 
give Lemma 3.10. □ 



4 Proof of Theorem 2.2 

4.1 Existence 

In this section, we assume that uq is a nonnegative function in L°°(Q) satisfying (5). For 

5 G (0,1), (us,(fs) denotes the solution to (KS)s constructed in Section 3. To prove 
existence of a weak solution, we use a compactness method. For that purpose, we first 
study the compactness properties of (us, fs)s- 

Lemma 4.1. There are functions u and ip and a sequence (b~ n )n>i, S n — > ; such that, for 
all T > and p G (1, oo), 

u 5n — ► u, in L p ((0,T) x U) as S n -> 0, (39) 

<p 5n — > <p, in Z7((0, T); W 2 ' P (Q)) as <5 n -> 0. (40) 
In addition, u G L°°((0, T) x 0) /or all T > and is nonnegative. 
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Proof. Thanks to Lemma 3.4 and Lemma 3.9, (it™) 5 is bounded in L 2 ((0, T);H l (f2)) while 
(d t uf)s is bounded in ^((O.T); (W 1 ' N+1 )' (ft)) by Lemma 3.10. 

Since H X (S}) is compactly embedded in L 2 (Sl) and L 2 (f2) is continuously embedded in 
(W^ N+1 )'(Q), it follows from [19, corollary 4] that (uf) is compact in L 2 ((0,T) x 0) for 
all T > 0. Since r i — )• is ^-Holder continuous, it is easy to check that the previous 
compactness property implies that (us) is compact in L 2m ((0, T) x Q) for all T > 0. There 
are thus a function n € L 2m ((0,T) x il) for all T > and a sequence (<5 n )n>i such that 

u Sn — ► n in L 2m ((0, T) x 0) as 5 n -> 0, (41) 

for all T > 0, owing to Lemma 3.9, we may also assume that 

u 5n A u in L°°((0, T) x 0) as 5 n -> 0. (42) 

for all T > 0. It readily follows from (41) and (42), and Holder inequality that (39) holds 
true. Since elliptic regularity ensure that 



V8 n \\w*>r ^ C (p) \\ u S k ~ u S n \\p, 



for all k > 1, n > 1, and p G (1, oo), a straightforward consequence of (39) is that (ifs n ) n >i 
is a Cauchy sequence in L p ((0, T); W 2,P (Q)) and thus converges to some function ip in that 
space. Finally, the nonnegativity of n follows easily from that of us n by (39). □ 



Proof of Theorem 2.2 (existence). It remains to identify the equations solved by the limit 
(it, ip) of (us n ,<P5 n ) n >i constructed in Lemma 4.1. To this end we first note that , owing 
to (39) and the boundedness of (ug n ) n and u in L°°((0,T) x fi), we have 



u m in L p ((0,T) x SI) as 5 n -> 0, 



(43) 



for all T > 0. Since (V(us n +S n ) 2 ) n > 1 and (Vm™ ) n >i are bounded in L 2 ((0, T) x fi) for 
all T > by Lemma 3.4 with p = 2 and p = m + 1, we may extract a further subsequence 
(not relabeled) such that 



m + 1 

v(u 5n + <y 2 



Vn— in L 2 ((0,T) x SI), 
Vn^ ->> Vn m in L 2 ((0,T) x 0), 
for all T > 0. Then if V £ i 4 ((0, T) x O; M w ), 



(44) 
(45) 







■T 



Vn m ] dxds 



< 



+- 



< 



+ 



m + 1 
2 

m + 1 
2 

m + 1 
2 

m + 1 
2 

m + 1 



V" 



T 



0<5n + S n 



m — 1 
1 2 



V(n 5n + <5 n ) 2 



m — 1 m+1 

u 2 Vn 2 



dxds 



771 — 1 , 

U 2 



dxds 



T 



m—1 

u 2 ip 



v(u 5n + 5 n y 



m+1 

Vn 2 



dxds 



4 l|V(n 5n + <5 n ) 



b ||(«5 n + <W 



m — 1 

u 2 U 



m—1 

u 2 ip 



V(li5„ + <5« 



7^^+l 



m + 1 

Vn 2 



dxds 



bmce n 2 



V> G L 2 ((0,T) x 0), we deduce from (39) and (44) that the right-hand side of 
the above inequality converges to zero as n — > 00. In other words, 



V(us n + S n 



Vu m in L*((0,T) x SI), 



(46) 
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for all T > 0. 



Now, we are going to show that (u, (p) in Lemma 4.1 is the desired weak solution in 
Theorem 2J. Let T > and i/j G C l {[®,T) x Q) with = 0. The solution of (7) 

satisfies 



f I [V(n 5n + 8 n ) m ■ Vip - u 5 „ V Wn ■ - u Sn d^} dxdt = [ 
Jo Jn Jn 



u Q tp(0,x) dx, (47) 



and, 

/ [ [V<^„ • Vip + M ip — u$ n ip] dxdt = 0. (48) 
Jo Jn 

From (46) we see that 

[ [ v (^5 n + $n) m ■ dxdt — > [ [ Vu m ■ VV> dxdt as 8 n ->■ 0. 
Jo Jn Jo Jn 

From (39) we get 

/ / u$ n dtip dxdt — > / u dtip dxdt as 5 n — > 0. 
Jo Jn Jo Jn 

From (39) and (40) we get 



/ ug n V(^,5 n • dxdt — > / u Vtp • Vip dxdt as 5 n — > 0. 
10 Jn Jo Jn 

Thus we conclude that u satisfies 

[ (Vu m ■ Vip - uVip ■ W - u ■ d t ip) dxdt = [ u (x) • ip(0, x) dx. 
10 Jn Jn 

Similarly, from (40) we see that 

/ Vip Sn ■ Vi/j dxdt — > I [ V<p-Vif> dxdt as 5 n 0, 
Jn Jo Jn 

and from (39) we see that 

/ u$ n ip dxdt — > / u tp dxdt as 5 n — > 0. 
Jn Jo Jn 

Thus, we have constructed a weak solution (u,(p) of (KS). □ 



4.2 Uniqueness 

In this section, we prove the uniqueness statement of Theorem 2.2 under the additionnal 
assumption (6) on ip. The proof relies on a classical duality technique, and on the method 
presented in [2] 

Proof. The proof estimates the difference of weak solutions in dual space H 1 ^)' of H 1 ^), 
motivated by the fact that the nonlinear diffusion is monotone in this norm. 
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Assume that we have two different weak solutions (ui,(p%) and (u2,<P2) to equations 
(1) corresponding to the same initial conditions, and fix T > 0. We put 

(u, <p) = («i - u 2 , tpi - (pi) in [0, T] x Q. 

Then ip is the strong solution of 



—Aip 
d u (p 
< ip > 



u in $7, 
on dQ., 
0. 



(49) 



Since d t u e L 2 ((0,T); H 1 ^)'), we have 

—Adtf = dtui — dtU2 = dtu in H 1 ^)', 

and 



ld_ 

2di 



liv^l 



V99 • Vdt<p dx 
= - < Ad t <p, (p >(my,m=< d t u, ip >( H ^y,m 
Now it follows from (1) that u satisfies the equation 



d t u 
d v u 
u(0, x) 

Substituting (51) in (50) . we obtain 
1 d 



div(V(«5" - uf)) - div(mVv3 + uV</> 2 ) 

0. 



(50) 



(51) 



2dt 



ivpii! 



/ (uf - uf) Aip dx + / ui |V^| 2 dx + u Vip 2 • V99 dx. (52) 
Jn Jn Jn 



The first integral on the right-hand side of (52) is nonnegative due to the fact that z 1— >• z m 
is an increasing function. The second integral on the right-hand side of (52) can be 
estimated by 

/ u\ |V</?| 2 dx < ||«i||oo / |V(/3[ 2 dx. 
Jn Jn 

For the last integral, using an integration by parts we obtain 
u Vip 2 ■ Vip dx = 



/ u Vip2 ■ V<p dx = — Aip V(f2 ■ V</2 dx 
Jn Jn 

= Vip- V(Vcp 2 • V<^) dx 
Jn 

= ^ I diV dfjip 2 djip dx + ^2 d if d jV2 d ij^ dx - ( 53 ) 

i,j ^ n i,j ^ n 

integrating by parts the second integral on the right-hand side of (53), 

y~] J dip dj(p 2 dfj(p dx = Yl\ J ^il^l 2 dx 

= I I Vv9 2 • V([Vv?| 2 ) dx 
2 Jn 

= -\j^V2 \Vip\ 2 dx 
< C(T)\\\/p\\l 



1 -j 
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since — A</?2 = U2— < ^2 L°°((0,T) x f2). Together with (53) the previous inequality 
implies 



/ it V(f2 ■ V(p dx 
Jn 



< C(T) [ (IDV2I + I) \V<p\ 2 dx. 

Jn 

< (ll^2||L°°((o,r);VK 2 .°°(m) + 1) / |V<p|" 

Jn 



dx, 



provided that the L°°((0, T); W 2,oc (Q)) norm of the function ip% is bounded. Thus, sub- 
stituting the above estimates in (52), one finally obtains 



d_ 

dl 



[ \Vip\ 2 dx < C(T) [ \Vip\ 2 dx. 
Jn Jn 



(54) 



Notice that [|Vy(0)[|2 = which follows from (49) and the property u(0) = 0. Thus, 
inequality (54) implies 

||V^)||Be C(T)i ||V^(0)||2 = 0. 

Consequently, S7(p(t) = for all t G [0, T] and, since < (p(t) >= 0, we have ip(t) = 
for all t £ [0, T]. Using (49), we conclude that u(t) = for all t € [0, T]. Consequently 

(Ul,(fl) = (U2,(P2). □ 
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